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© A process and related apparatus for producing at least two matched color displays of a digital image using 
two different display devices (12,14), comprising using an adaptor (22) to convert the digital information 
representing the image to digital information such that the displayed image as a result of this converted digital 
information on one of the devices, appears the same as the image displayed on the other to an average 
observer when viewed under similar viewing conditions. The adaptor contains a transform LUT (28) constructed 
using display device modeling to generate corresponding values to preselected input digital image values to 
map the input color values Into the color space of the display device so that image matching is achieved. 
One or both display devices (12,14) can be replaced by one or two printers. 
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COPYRIGHT NOTICE 

A portion of the disclosure of this patent document contains nnaterial which is subject to copyright 
protection. The copyright owner has no objection to the facsinnile reproduction by anyone of the patent 
5 document or the the patent disclosure, as it appears in the Patent and Trademark Office patent file or 
records, but otherwise reserves all copyright rights whatsoever. 

BACKGROUND OF THE INVENTION 

10 1. Field of the Invention. 

This invention relates to color matching and more particularly to a process and related apparatus for 
matching the color displayed on a plurality of color display devices. 

75 2. Description of Related Art. 

The primary use of this invention is in the field of computer assisted color publishing systems 
particularly in the area of color matching. In such systems typically a color image Is scanned using a 
scanning device which measures light intensity reflected or transmitted from the image to generate 

20 electronic signals representing three achromatic images of the original, each of the three achromatic images 
resulting from detecting the light intensity from the original color image through a colored filter. The filters 
used are almost always Red. Green and Blue. The electronic signal from the scanner is converted to a 
digital signal in which light intensity levels are represented as numbers. Information identifying each set of 
numbers representing the image information obtained through each filter is also preserved. Thus, through 

25 these steps, the original colored image is converted to a plurality of image values, and for each picture 
element in the picture, there are three such values, a Red, a Green and a Blue. 

The image represented by the image values may be displayed in a CRT type monitor, or may be 
printed using a printing device able to accept electronic input. Often the image will be displayed at different 
locations and times using more than one display device, such as a plurality of CRT displays for observation 

30 and study by more than one operator. Similarly, hard copies may be desired in two or more different 
locations using two or more different printers. Even though the input values to the multiple CRT displays are 
the same, the displayed colors are vastly different, as anyone is well aware of who has ever observed the 
multiple television displays in a typical television sales store front. However, when one needs to make 
decisions about acceptability of color for a display, one needs to know with a great degree of confidence 

35 that the color, that one Is observing and discussing with an operator who is observing the same image on a 
different CRT, is the same as the color observed by the operator. The same is true, if the displayed image 
is one created on a printer and compared with the same image created on a different printer. 

Colorimetry, which is the study of color based on both spectral distribution of the energy reflected or 
transmitted from a test sample and the response of the human eye. as well as the spectral distribution of 

40 the illuminating source, provides a method to describe and measure color and enables one to determine 
when colors match. Through the use of CIE defined Tristimulus Values (which are the amounts of three 
primary lights which when added produce a visual, or colorimetric match with an original color), one may 
determine with reasonable certainty that if two colors have the same three CIE Tristimulus values that is. if 
the Red Tristimulus value of one color is the same as the Red Tristimulus value of the other, and so on for 

45 the Green and Blue Tristimulus values, then the appearance of the two colors will be the same to the 
average observer. CIE stands for the International Commission on Illumination. 

It appears, therefore, that the problem of matching color outputs of different displays is readily resolved 
by providing displays having the same CIE tristimulus values for the same Image value inputs. 

While the solution in principle is simple, creating a conversion device which will make two or more 

50 different displays produce the same tristimulus output values for the same input image values is very 
difficult. Each of the display devices operates in its own, device dependent, color space where image 
values at its input are transformed into display image values. The difficulty lies primarily in matching the two 
transformations occurring within the two display devices for the image values appearing at the input of 
each, so that both display the same displayed colorimetric values for the same input image values. 

65 The prior art solutions to color matching fall into two fundamentally distinct approaches. The first is 
based on the decomposition of a color vector to a set of primaries, and known as the primary decomposi- 
tion technique. Typical of this approach is the technique described in U.S. Patent application 07/494.463 
filed March 16. 1990. assigned to E. I. du Pont de Nemours and Company. According to this technique, the 
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input Image is decomposed into the linear combination of eight primaries (yellow, magenta, cyan, black, 
red, green, blue, and three-color black). Coefficients for a 4x8 matrix (the values of CYMK for each of the 
eight primaries) are adjusted such that the color of the eight primaries in the input from both systems to be 
matched are matched, and the new input values are found by matrix multiplication of an 8X1 matrix 

5 (coefficient of decomposition) with the 4X8 matrix. This technique suffers because of the non linearity in the 
additive properties of color dyes and because the "primaries" used are not mathematically Independent. 

The more common solution, adopted by the printing industry is the grid sampling technique. This 
involves using a color transformation formula such that the error between the targeted color and the 
processed color Is minimized. The differences between the various methods in existence are found in the 

10 specific transformation formulas. U.S. Patent No. 4,500,919 issued to Schreiber is a good example of the 
use of a transformation formula. The difficulty with this approach again lies in the non-linearity of the color 
addition process and the complexity of the color surfaces in a set of equations. 

SUMMARY OF THE INVENTION 

76 

The present invention relates to provide a method for matching the color display of at least a first and a 
second display devices, comprising: 

(I) creating a transform LUT for converting Input color values to output color values by: 
(1) producing a first preselected plurality of input color values; 

20 (2) using said plurality of input values to display a plurality of color patches in said first and said 

second display devices; 

(3) obtaining a colorimetric value of each of the displayed color patches in each of the displayed 
devices and using said colorimetric values to create a first and a second model for the first and 
second devices respectively correlating preselected input color values to displayed colorimetric values 

25 for each of the two devices; 

(4) inputting to said first and second models a second preselected plurality of color values; 

(5) comparing the output of the first model to the output of the second model to obtain an error signal 
indicative of the difference between the two output signals; 

(6) using the error signal to modify the Input color values to the second model and again comparing 
30 the output of the first model to the output of the second model to obtain a new error signal; 

(7) repeating the process of steps (5) and (6) above until the error signal is a minimum; and 

(8) using the modified color values to create a transform LUT correlating input values to modified 
values; 

(II) using the transform LUT to transform any plurality of source color image values before inputting said 
35 image values to an Input of the second display device; and 

(III) displaying said transformed source color image values on said second display device. 

When there is no corresponding source Image value In the transform LUT correlating an input source 
color image value to a modified color image value, one can use interpolation to derive a corrected source 
image value from a closest source and corresponding modified value in the LUT. 
40 Preferably, the first and second preselected pluralities of color Image values are the same. 

This invention further relates to an apparatus for generating a transform LUT for converting input color 
values to output color values comprising: 

(1) digital color image values input means; 

(2) first and second display model means for producing each an output colorimetric tristlmulus Image 
45 values for input color values, each of said model means having an input and an output, the second 

model input connected to the input means; 

(3) an adder device having a first signal input connected to the Input means, a second, correction error 
signal input, a control signal input, and an output connected to the first model input, for outputting 
modified color values; 

60 (4) means connected to the output of said first and second model means, for comparing the tristlmulus 
values output of the first and second display models and for producing an error signal; 
(5) means for testing the error signal to determine if said signal is a minimum error signal, and for 
outputting a correction error signal and a control signal to the adder, said testing means connected 
between said means to compare and said adder device; and 

55 (6) means also connected to the Input means, the adder output, and the means for testing, for receiving 
the adder output and for generating a transform LUT correlating digital input color values to the modified 
color values appearing at the output of the adder when the means for testing determines that the error 
signal is a minimum. 
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The apparatus may further comprise means for storing the generated LUT. A gamut mapping means 
may be included between the output of the second model and the means to compare for mapping output 
values from the second model into a color gamut commensurate to that of the first model. 

The apparatus may comprise hardware, or may be a computer programmed through software to 
6 perform all of the above operations, or may be a combination of dedicated hardware and computer 
implemented software. 

BRIEF DESCRIPTION OF THE DRAWINGS 

70 The invention can be more fully understood from the following description thereof in connection with the 
accompanying Rgures described as follows. 

Figure 1 is a schematic representation of an arrangement In which two matched images are displayed 
on two CRT displays. 

Figure 2 is a schematic representation showing the arrangement used in obtaining data for generating 
75 the LUTs used in modeling the two displays of Figure 1. 

Figure 3 is a schematic representation of the process for obtaining the transform LUT used in the 
adaptor in the arrangement shown In Figure 1 . 

Figure 4 is a generic representation of a display model. 

Figure 5 is a schematic representation of the process used in deriving the correction factors shown In 
20 Figure 3. 

Figure 6 Is a schematic representation of an arrangement in accordance with this invention in which two 
matched images are printed using two printing devices. 

Figure 7 is a schematic representation showing the arrangement used in obtaining data for generating 
the LUTs used in modeling the two printing devices of Figure 6. 
25 Figure 8 Is a schematic representation of the process for obtaining the transform LUT used in the 
adaptor In the arrangement shown In Figure 6. 

Figure 9 is a schematic representation of an arrangement in accordance with this invention in which two 
matched images are printed using a CRT display and a printer having including an RGB to CMYK 
converter. 

30 Figure 10 Is a schematic representation of an arrangement in accordance with this invention in which 
two matched images are printed using a CRT display and a color proofing device using a YMCK to RGB to 
drive the CRT display. 

Figure 1 1 is a generic representation of a display model for use In developing the transfer LUTs 
required for the embodiments shown in Figures 9 and 10. 

35 

DESCRIPTION OF THE PREFERRED EMBQDIMENT(S) 

The invention will next be described in detail with reference to the drawings in which similar characters 
indicate similar elements in all figures of the drawings. During the following discussion, color values are 
40 quantized in 256 steps corresponding to an 8 bit system. Other quantizations are possible, and not intended 
-> to be excluded because of the use of the 256 steps in the following examples. 

The invention comprises a method and apparatus for providing color visual matching of two representa- 
tions of the same image when the image is displayed in two different displays, so that an observer will on 
visual examination of the displayed Images observe two Images having substantially the same color 
45 appearance. Using the method or apparatus of this invention causes each of two observers, whose visual 
response substantially conforms to the visual response of an average observer as defined by the CIE 
Institute, upon observing any one of the color displays under the same or similar surrounding illumination 
and background, to receive substantially the same visual impression. 

Figure 1 shows in schematic representation a situation where a color image is displayed on two CRT 
50 type displays which receive RGB type inputs to display an image. The Image may be stored and 
manipulated in digital format in a work station 10 which itself may comprise a scanner, an image processor, 
a display device, and have image storage capabilities, or any portion of the above, or more. 

Because the required input to displays 12 and 14 in this example is a digital RGB type signal, the work 
station outputs over line 16 a digital RGB signal representing a color image. The digital RGB signal 
55 comprises a set of color values, one set of three values (R, G, B) for each picture element of the image to 
be displayed. The signal is directed over line 18 to the first display monitor 12 on which the colored image 
Is displayed. The digital RGB signal is also sent over line 20 to the second display monitor 14, which could 
be located at a location different than the location of the first display 12. 

4 
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Prior to applying the RGB signal to the second monitor 14, there Is interposed an adaptor 22 which 
receives the digital RGB signal and after processing the signal, outputs a signal comprising a new set of 
digital color values R',G',B* over line 24 to an input of the second display monitor 14 for each image pixel. 

The adaptor 22 comprises a receiving and determining element 26 which receives each set of R,G,B 
s values and checks them against a plurality of sets of R,G,B values comprising a transform look up table 
(LUT) 28 which correlates R,G,B values to a new set of R',G*,B' values which when used as an Input to the 
second display monitor 14 will generate a pixel whose colorimetric value will match that of the pixel 
generated by the same R,G,B set of color values In the first display monitor 12. 

If the receiving and determining element 26 identifies a corresponding set of R.G.B values in the LUT 
10 28. it outputs the appropriate R\G\B* values on line 24. If not. the R,G,B input is directed to an interpolator 
30 . where R',G\B* values are derived by interpolation using existing adjacent sets of R.G.B values in the 
LUT. Preferably, the selection of the R,G,B values for the development of the transform LUT 28 is such that 
values derived by Interpolation produce a visually acceptable match of the displayed pixel in the two 
monitors 12 and 14. 

75 For this system to operate with any degree of success, the development of the transform LUT 28 is 
critical. The values generated by this LUT 28 must indeed lead to values that correctly reproduce colors on 
the second display 14 which visually match the colors on the first display 12. Figures 2. 3 and 4 help 
explain how the LUT 28 is generated. 

The first step In developing the transform LUT 28 is the creation of two models representing the two 

20 displays whose output is to be matched. Figure 2 shows how this is accomplished. 

The work station 10 is used to generate and output over line 16 a set of digital R.G.B color values in a 
regular pattern which is selected to supply color values closely enough spaced so that interpolation of 
values between existing values is reasonably accurate. Preferably, but not necessarily, the output is 
adjusted so that the color values produce displays 32 and 34 on the first and second display monitors 1 2 

25 and 14, respectively, comprising a plurality of different color patches on display screens of the display 
monitors 12 and 14. As an example, the work station 10 may generate a plurality of digital RGB values such 
as shown in the following table I: 

TABLE I 

30 



R: 


0, 


13, 


26, 


51, 


76, 


102, 


128, 


153. 


178, 


204. 


230. 


255. 


G: 


0, 


13, 


26, 


51, 


76, 


102, 


128. 


153. 


178. 


204, 


230, 


255. 


B: 


0. 


13. 


26. 


51. 


76. 


102. 


128. 


153. 


178, 


204, 


230. 


255. 



35 

Each combination of R.G.B values from this table represents a set (R,G,B)n of color values. In this 

illustration where there are 1728 possible such sets of R,G,B. values, n = 1 to 1728. 

These 1728 sets of R.G.B. values are supplied to both display monitors 12 and 14 over lines 18 and 20' 

creating a total of 1728 displayed patches for each of the monitors 12 and 14, corresponding to the 1728 
40 different combinations of the selected values. All 1728 patches are not necessarily displayed on each 

monitor screen simultaneously. 

A colorimetric measuring device 36 Is used to read the displayed patches displayed in the first and 

second monitors 12 and 14 as shown In Figure 2. The output of the colorimetric device is a colorimetric set 

of color values for each of the patches. In our preferred embodiment, the output of the colorimeter are CIE 
45 defined L,a and b sets of color values in the Lab color space. Output colorimetric measurements given in 

other colorimetric color spaces are equally acceptable, such as XYZ, Tristimulus CIE defined RGB, etc. 

Uniform as well as non-uniform color spaces may be used. However, values in non-uniform color spaces 

related to a uniform color space with a known mathematical relation may result in increased calculation 

steps which must be compared with any gained advantage as a result of this selection, to determine if the 
50 particular choice of color space is justified. (For a description of the different color spaces and related 

terminology, see. in general. DIGITAL IMAGE PROCESSING. 2nd edition, by William Pratt, published by 

John Wiley and Sons, Inc. pages 62-73.) 

Once the measurement of all patches is completed, two LUTs are compiled, representing the transfer 

functions of the two monitors 12 and 14. The first LUT will consist of the (RGB)^ sets of values and the 
55 corresponding (Lab)n sets of values read off the display on monitor 12 and the second LUT will consist of 

the same (RGB)n sets of values and the corresponding (L'a'b*)n sets of values read off the display monitor 

14. 



5 



>JS0OCID: <EP 



0604755A1 I > 



EP 0 604 755 A1 



Referring to Figure 3. a first and a second model 40 and 42 are built representing the two monitors 12 
and 14, respectively. Figure 4 shows a generic model 44 structure, which is used throughout this 
description in this invention to convert color values. The model 44 comprises an LUT 48, which Is the 
particular device derived LUT representing the device transfer function. In this example, the LUT 44 for the 
5 first model 40 will be the first LUT correlating the (RGB)n sets of color values to the (Lab)n sets of values, 
and the LUT 44 for the second model 42 will be the second LUT correlating the same (RGB)n sets of color 
values to the (L'a*b')n sets of values. 

In addition to the LUT 48, the model 44 comprises a receiving and determining means 46 similar to, or 
the same as. the receiving and determining means 26 described earlier, and a mathematical interpolating 
70 means 50, also similar to, or the same as, the previously described interpolating means 30. 

Figure 3 schematically represents the generation of the transform LUT. Sets of R,G,B color values from 
a preselected plurality of R.G.B sets of color values are used. This plurality of R.G.B sets of values is. 
preferably, also produced in the work station 10, and, again preferably. Is the same as the sets of (R,G.B)n 
values from Table I previously used to generate the two LUTs for the two display models 40 and 42. 

75 Each (R.G,B)n set of values is directed to the input of the first model 40 over lines 54 and 58. Model 40 
produces an output of (Lab)n values corresponding to the Input (R,G,B)n values for this model 40. The same 
(R.G,B)n values are directed to the second monitor model 42 over line 56. Ahead of model 42. there is an 
adder 60 which operates to add to the Rn. Gn. and Bn components of the (R,G,B)n set of values any 
correction factor dRn. dGn. dBn, appearing on line 78. At first, nothing appears on line 78, and the Rn, Gn 

20 and Bn values are applied to the input of model 42 unaltered. Model 42 also produces an output set of 
(L'a'b')n color values corresponding to the input (R,G,B)n. This output appears on line 64. 

The (Lab)n values over line 72 and the (L*a'b*)n values over line 64 are next compared in comparator 66 
and a difference signal [(Lab)n-(L*a'b*)n] is generated and directed over line 74 to correction factor generator 
76. Using this difference signal, correction factor generator 76 produces correction factors dRn(i), dGn(i), 

25 dBn(i) (in a manner to be explained herein below) which are sent over line 78 and added to the 
corresponding components of the (R,G,B)n set of values in adder 60, resulting in a new set of color values 
Rn(i)* = Rn + clRn(i), Gn(i)' = Gn + dGn(i ), and Bn(i / == Bn + dBn(i ) on line 62. These new values are again applied 
to the input of model 42 which produces a new output set of L2'a2'b2' values on line 64 which is again 
compared in comparator 66 with the set of (Lab)n values from model 40. A new set of correction factors 

30 dR„(2,, dGn(2). and dBn(2) is produced and added to Rn(i)', Gn(i/ and Bnd/ to produce a new set of Rp'. G„', 
and Bn* values. An error "E" defined as: 

E = [(L-U)2 + (a-a*)2 + (b-b*)2]^ 

35 is used to determine when to stop this cycle. Referring to Figure 5. error "E" is compared In comparator 92 
with preselected minimum acceptable limits or a "0" value. If the error is "0", or within the preselected 
acceptable limits, a signal over line 94 switches switch 98 (which is part of adder circuit 60) to feed the 
Rn'.Gn.'and Bn* values which produces the minimum over line 80 to the transform LUT compiler 82, rather 
than to model 42; the next input set of values (R,G,B)n+i is then applied over line 54 to the two models and 

40 the whole process is repeated for this new set of values and so on. 

The set of Rn'Gn'Bn* values which produced this zero or acceptable minimum error are sent over line 80 
to the LUT compiler 82 where the (R.G,B)n and corresponding (R*G'B')n sets of values are used to compile 
the transform LUT 28 so that for each set of (R,G.B)„ values there is a corresponding set of (R',G'.B*.)n 
values. Means to store the transform LUT 28 is Included in the compiler 82. This LUT 28 is the the same 

45 transform LUT shown in Figure 1 discussed earlier. 

The correction factor generator 76 produces the dR, dG and dB correction factors as follows. The input 
to the correction factor generator 76 is the difference signal produced by the comparator 66, specifically [- 
(Lab)n-(L'a*b')n]. The individual components of this signal are: dL which equals (L-L'), da which equals (a-a') 
and db which equals (b-b'). The correction factors dR, dG and dB, and the difference signals dL, da. and db 

50 are related by the following relationships: 

(1) (db/dG)da-(da/dG)db 
= [(db/6G)*(da/dR) - (da/dG)*(db/dR)]dR 
+ [(db/dG)*(da/dB) - (da/dG)*(db/dB)]dB 
55 (2) (da/dG)dL-(dU6G)da 

= [(da/dG)*(dL/dR) - (dL/dG)*(da/dR)]dR 
+ [(da/dG)*(dUdB) - (dLydG)*(da/dB)]dB 
(3) dL = (dL/dG)dG + (dL/dR)dR + (dL/dB)dB. 

6 
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The actual numeric values for dUdR, 6LydG, dlVdB, etc. for use in solving the above system of 
equations, are obtained by the process schematically illustrated in Figure 5. The input Rp.Gn.Bn values are 
changed by holding two the same and incrementing the third by one unit in incrementing module 96. This 

5 module comprises a memory for temporarily holding the Input values of Rn, Gn, and Bp and an arithmetic 
means for incrementing each of those values by 1 and for applying various combinations of Rn. Gn. Bn, 
Rn + 1 . Gn + 1 , and Bn + 1 sets of values to the input of model 42. 

An example which uses illustrative RGB values, will be used to explain the derivation of the numerical 
values for the partial derivatives needed to solve the above equations. (1) (2) and (3). Let the input Rn, Gn. 

10 and Bn values for a set of (R,G,B)n values be 100, 100. 100: Assume that when these values are applied to 
both models 40 and 42, two sets of Lab values, (Lab)n and (L'a'b')n are produced such that the error E as 
previously defined is different than "0" or a preselected acceptable minimum. Incrementing module 96 
produces over line 65 three new color value sets. Rn,Gn,(Bn + 1). having the numerical values (100. 100, 
101) in this illustrative example. Rn,(Gn + 1).Bn, having the numerical values (100.101,100) and (Rn + 1)- 

75 .Gn.Bn, having the numerical values (101.100,100). These new values are applied to model 42 over line 63'. 
The resulting variations in the components of the (L*a'b')n values produced by model 42 are directed to an 
arithmetic calculator 82 over line 75 and to buffer memory 84. The calculator 82 performs the operations: 
(L-L')/R-(R + 1), (L-L')/G-(G + 1), (L-U)/B-(B + 1). (a-a')/R-(R+ 1). (a-a')/G-(G + 1). (a-a*)/B-(B + 1 ). (b-b')/R- 
(R + 1), b-bVG-(G + 1) and b"bVB-(B + 1), to derive the partial derivative values (dL/dR. dUdG. etc.) used in 

20 equations (1), (2) and (3) above. 

The results of the calculator 82 are supplied over line 83 to another calculating unit 86. The dL. da. and 
db values from comparator 66 are also directed to this calculating unit 86. Unit 86 solves the system of 
equations (1), (2) and (3) for dR, dB and dG to derive the correction factors. 

A problem Is sometimes encountered in determining a set of values that will produce a best color 

25 match on two different displays for the same input set of color values. The problem is that a minimum error 
may be arrived at, which may not be an absolute minimum, but only what may be referred to as a local 
minimum, and there may be other combinations of values (R', G\ 3% which could produce a smaller even 
a zero error and a better color match. To test for this possibility, an optimization operation is performed. A 
buffer memory 90 holds the input (R,G,B)n values, the corresponding (R\G',B*)n values and the error "En" 

30 resulting from the use of this combination. The next set of (R.G.B)n+i values is processed and it, the 
corresponding (R'.G',B.)n+i values, and resulting error "En", are also stored. The following (R.G.B)n+2, 
(R'.G',B')n+2 and "E3" are stored, etc., until a set of (R,G,B)n+x and (R'G'B')n+x values for which the error 
"Ex" = 0 Is encountered. Once this true "0" error has been identified, a backtracking operation is performed 
to determine if a "0" error can be obtained by correcting the set of (R,G,B)n«x-i values using as an input to 

35 the model 42 the (R.G,B)n.x set of values (which is the set for which a "0" error was obtained) while model 
40 receives the (R,G,B)n-x-i set. This operation will produce a new set (RI'.GV. B1*)n-x-i of values and a 
new error "E1Vx-i for the input values (R.G,B)n_x-i. If "E1"n-x-i = "0", the (Rr,Gr,Br)n.x-i set of 
values is substituted for the (R'G'B')n-x-i set of values, and the transform LUT 28 is recompiled using these 
(Rr,Gr,B1')n=x-i values as corresponding to the (R.G.B)„.x-i seL 

40 If no "0" error Is obtained with the above process, then the earlier arrived at result minimum is 
compared and if lower than the newly found minimum, the result stays. Otherwise the R\G\B' values giving 
the least error are used in compiling the transform LUT 28 as they will result in a more accurate match. The 
optimization process Is then repeated for the next (moving backwards) set of (R.G.B)n-x-2 input values and 
sequentially for all (R,G,B)n sets of values for n<x always using the values from the set with "Ey" = 0, 

45 whose ny number Is the closest to the tested set n number. When all stored sets of (R,G.B)n values have 
been tested in this manner, the memory contents of the buffer memory 90 are erased and a new set of 
input (R,G.B)n beginning with the next n, i.e.. n=x + 1 is applied over line 54 and the process of compiling 
the transform LUT 28 continues. 

For Illustration purposes, lets assume that an input set of values (R,G,B)2o = (100, 100, 100) is applied 
50 to both the first model 40 and the second model 42 and that this results In a corresponding set of values 
(R',G'.B')2o = (100. 102, 103,) which gives a minimum acceptable error but different from zero. Let us 
further assume that the next set of input values Is (R,G,B)2i and that (R,G,B)2i = (110, 100, 100). Also let 
us assume that the corresponding set of values (R'.G*.B')2i are (105, 100, 101), and that the error in this 
instance is zero. 

65 When the zero error is detected, the optimization process begins, and no new input sets of (RGB)n 
valuesare applied to line 54. 

The set preceding the set which had the zero error, i.e., set of input values (R,G,B)2o = (100, 100, 100) 
is applied to the input of model 40. However, the set (R\G',B')2i = (105. 100, 101) Is applied to the Input of 
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the second model 42 (rather than the set (R,G.B)2o) and a new set of corresponding (R1*G1'B1')20 values is 
arrived at, i.e., 101, 102, 103. which will either give a zero error or another minimunn error. If the result is a 
zero error or jf the new error is smaller than the previous, the new (RVGVBV)2o values are substituted for 
the original values (R'G'B')2o that were used earlier in the compilation of transform LUT 28 and will produce 
5 a closest match. 

The process described is best performed in a computer, or if in hardware, under the control of a 
computer suitably programmed to perform the above computations. The diagram of Figure 5 is provided for 
illustration, rather than as an actual working embodiment of hardware sufficient to perform the above 
process, so as to enable one to properly program a computer for this task. 

10 Referring again to Figure 3, there is shown a block 70 between the output of model 40 and the input to 
comparator 66. This represents an optional color gamut correcting device, which may be used to match the 
output color gamuts of the two display devices. Gamut matching is often often used because the color 
range that is possible to reproduce with a particular display is usually not identical with the color range 
possible to reproduce with another. Thus, while both displays have a large number of colors that are 

75 reproducible by both, there may be certain combinations of sets of color values that are outside the 
reproduction capability of one or the other display. While the use of the transform LUT 28 in itself inherently 
provides color gamut matching between the two devices considered, additional color gamut matching may 
be used In the development of the LUT 26, by inserting a gamut matching device 70 as shown. This device 
will perform a preselected mathematical manipulation to the (Lab)n set of values so that there will not be 

20 any sets for which there are no corresponding (L'a'b')n values from the model 42. Such mathematical 
manipulation may be as simple as simple truncation or as complicated as data compression. 

These types of mathematical operations are well known in the art of color matching. For a brief 
discussion of gamut matching and the necessity of color compression, also see "COLORIMETRY FOR 
ELECTRONIC IMAGING SYSTEMS, An intensive short course for Scientists and Engineers", Copyright 

25 1991, Munsell color science laboratory, Rochester Institute of Technology, by R. S. Berns and R. S. Hunter, 
Chapters 27 and 28. 

Figures 6, 7 and 8 illustrate the case where the images that are to be matched are images produced in 
color proofers or printers which use a typical four color reproduction system with four color inks, such as 
Yellow (Y), Magenta (M), Cyan (C) and Black (K). 

30 Referring to Figure 6. as in the previous example, a work station 10 generates digital color values 
representing pixels of a color image. Since a four color system is used and the reproducing equipment 
accepts YMCK signals, the work station 10 output is a set of YMCK signals over output line 116. The 
signals are directed to a first printer 112 over line 118 and to a second printer 114 over line 120. Ahead of 
the input to printer 114 is placed an adaptor 122 which again comprises a transform LUT 128, a receiving 

35 and determining element 126 and an interpolator 130. The LUT 128 is a transform LUT which correlates a 
plurality of YMCK values to a second plurality of Y'M'C'K* values such that colored pixels produced using a 
set of YMCK values in printer 112 and a corresponding set of Y'M'C'K' values from the LUT 128 in printer 
114 will be identical or virtually identical to the standard CIE average observer when viewed under similar 
illumination and surroundings. 

40 As for the previous example the interpolator 130 is used to supply corresponding values for YMCK input 
values not listed In the transform LUT 128. The transform LUT 128 has been compiled using substantially 
the same process as for the transform LUT 28 used in the previous example. As illustrated in Figure 7, the 
work station 10 may generate a plurality of digital YMCK values such as shown in the following Table H: 

45 TABLE II 



Y: 


0. 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204, 


230, 


255. 


M: 


0. 


13, 


26, 


51, 


76, 


102, 


128, 


153, 


178, 


204. 


230. 


255. 


C: 


0. 


13, 


26, 


51, 


76, 


102, 


128. 


153, 


178. 


204, 


230, 


255. 


K: 


0, 


13, 


26, 


51. 


89. 


128. 


178. 


255. 











Each combination of Y.M.C.K, values from this Table II represents a set (Y,M,C,K)n of color values. In 
this illustration, there are 13824 possible such sets of Y,M,C,K, and n = 1 to 13824 values. 
55 The 13,824 patches resulting from the different color value combinations above are printed by both 

printers 112 and 114 to provide two sets 132 and 134 of 13.824 patches each. Again as before, a 
colorimetric measuring device 36 is used to read all of the above patches 132 and 134 and the colorimeter 
output is used to compile two LUTs representing the transfer functions of the two printers 112 and 114. 
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Again it is preferred that the colorimeter output is In the Lab color space and the LUTs correlate (Y,M,C,K)n 
inputs to (L,a,b)n outputs for- each of the printers 112 and 114. These LUTs are used to construct models 
140 and 142 using the generic model diagram shown in Figure 4 for the two printers 112 and 114. 
respectively. 

5 Figure 8 Illustrates the process followed for the compilation of the transform LUT 128. This process is 

similar to the process described earlier for the compilation of the transform LUT 28. The LUT, however, 
correlates two pluralities of color sets, each set having four rather than three variables. Of the four 
Cn.Yn.Mn.Kn Variables, the three Cn.Yn.Mp variables are independent variables, but the fourth, Kn, Is not. Kn Is 
dependent on the other three from which it is derived. The relationship between CYM and K is well known 

10 in the graphic arts. 

Preselected sets of (C,Y,M,K)n values from the work station 10 or elsewhere, are directed to the inputs 
of both model 140 and 142 over lines 154. 162, 156 and 158. It is prefered that the preselected values be 
the same as the values used to compile the LUTs used in the models 140 and 142 and for purposes of this 
description the same value combinations resulting using the values in Table 11. 

75 Each input set of values (C,Y,M.K)n is applied directly to the input of model 140 over lines 154 and 158. 
The same set of values is also applied over line 156 to one input of adder 160. The adder output is directed 
to switching circuit 198. The output of switching circuit 198 Is controlled by a signal over line 194 and may 
be directed to the Input of model 142 or to transform LUT compiler 182 over lines 162 and 180, 
respectively. The input set of values is also applied to compiler 182 over line 155. 

20 The adder 160 has a second input on which are applied correction factors dCn, dYn and dMn produced 
in the correction factor generator 176. When the (Y,M,C,K)n set of values first appears at the adder 160 
input, there are no correction factors applied to the adder and the values Yn.Mp.Cn.Kn are applied to the 
input of model 142 unaltered. Model 140 produces an output set of values (L,a,b)n and model 142 produces 
another set (L',a',b')n. These two sets are compared in comparator 166. For this description, it is assumed 

25 that an optional gamut adjustment module 170 is not used to modify the {L,a,b)n output of model 140 prior 
to applying It to the comparator 166 over line 172 

The comparator 166 produces a difference signal [(Lab)n-{L'a'b')n]. This output is directed over line 174 
to correction factor generator 176. Using this difference signal, correction factor generator 176 produces 
three correction factors dYn(i), dMn^ij, dCn(i) (in a manner similar to the manner in which correction factor 

30 generator 76 produced the dR, dG, and dB factors as previously explained). The correction factors are sent 
over line 178 and added to the corresponding components of the (Y,M,C,K)n set of values in adder 60. 
resulting in a new set of color values Y'n(i) = Y„+dY„(i,. M'„<i) = M„ + dM„<i,. and C'n(i) = Cn + dC„(i) on line 
162. 

The fourth value, Kn is stored in element 177 which comprises a buffer memory and a means to modify 
35 the stored value Kn by predetermined increments. At this time Kn remains constant and these new Y'n(i), 
M*r»(i), C*n(i), and Kn values are again applied to the Input of model 142 which produces a new output set of 
L2'a2'b2' values on line 164 which is again compared in comparator 166 with the set of (Lab)n values from 
model 140. A new set of correction factors dYn(2). dMn(2). and dCn(2) is produced and added to Y'no). M'no) 
and C'n(i) to produce a new set of Yn'. Mn*. and Cn* values. Again, an error "E" Is defined as: 

40 

E = [(L-L')2 + (a-a')2 + (b-b*)2]^'2 

The correction factor generator 176 also includes means to calculate, compare and store this error E„ 
and the corresponding (YMCK)n and (rM'C'K)n set of values that resulted in this error 

45 The above process is repeated continuously correcting the original input values Yn.Mn.Cn. with K^ held 
constant until the error E is zero or a predetermined acceptable minimum. If a zero error is found, then the 
two sets of values produce a perfect match and the search terminates. Switch 198 is actuated and the two 
sets, of values are used in the compiler 182 to compile the transform LUT 128. If only a minimum is arrived 
at It is again stored In the buffer memory and Kn Is incremented in unit 177 and the process repeated in 

50 search of another set of (Y"M"C"K')n values yielding another minimum. This new minimum Is compared 
with the previously found one. The (Y*M*C'K)n values or the values (Y"M"C"K*)n which produced the best 
minimum are kept. 

The correction factor generator In this instance solves a system of differential equations that is of the 
same format as previously discussed. The equations only differ in terms of variable designation and are as 
55 follpws: 

(1 a) (db/dY)da-(da/6 Y)db 

= [(6b/dY)*(da/dM) - (da/dY)*(db/6M)]dM 
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+ [(db/dY)*(da/dC) - (da/6Yndb/dC)]dC 
(2a) (da/dY)dL-(6L/dY)da 
= l(da/bYy{dUdWi) - (dL/dY)*(da/dM)]dM 
+ [(da/6YndL/dC) - {dUdCy(da/6C)]6C 
5 (3a) dL = (dLydY)dY + (6Ly6M)dM + (6UdC)dC. 

The partial derivatives are calculated in the same manner as the partial derivatives were calculated for 
the R.G.B. case, by Incrementing the Yn»Mn,Cn values by one, as explained earlier with the help of Figure 5. 
The issue of whether a minimum error En is a true minimum error or a local minimum, meaning that 

10 there may be another smaller error resulting from a different combination of (C',Y',M',K')n values which 
could be found by starting with a different input set to the second model 142, rather than the (Y,M,C.K)n set 
that was used is again a problem. The same testing and optimization process described earlier can be used 
In this case too, to test and discover whether there may be other better combinations of (YMVI'.C\K*)„ 
values, by holding in a buffer memory all the consecutive sets of (Y,M,C,K)n. the corresponding (Y'M*C'K)n 

75 sets and error £„ as described earlier. However, it has been found that more accurate results may be 
obtained by following a somewhat more calculation intensive process. This process utilizes two input sets of 
values which have resulted in a zero error as inputs to obtain two optimized sets of (Y1MVI1',C1\Kr)n and 
(Y2MV12',C2M<2')n values in addition to the (r,M*,C',K')n set and then selecting the set that produces the 
best match, i.e., least error, from the three sets. 

20 Because of the presence of the fourth color K, which is not an independent color as the Y.M,C, and 
which by its dependency on the other three can introduce a number of color combinations resulting in 
matching colors, the existence of two different sets of values both producing a perfect match is also 
possible. Further, the optimization process must be able to rationally select the set that will produce the 
better overall match if used in the compilation of the transform LUT 128. When two sets produce perfect 

25 matches, then the one that produces better interpolation results for input values that are not present in the 
LUT will be selected. 

The optimization process adopted for this case. Is based on selecting a predetermined fixed number 
"x" of sets of values (Y,M,C,K)n, the corresponding (y\C\M\K% and the associated error Epx and, then, 
proceeding with the optimization calculation using zero error resulting set values from either side of a target 

30 set to derive two sets of possibly lesser error resulting values in addition to the set originally found. Thus, 
the Input set (Y,M.C.K)n8, for Instance, for which the set of (Y'C'M'K')n8 gave us originally a minimum but 
non zero error, will be optimized using preceding and following values, to obtain (YV. CI', M1',K1')n8 values 
with data from a set (Y',CMVI',K%y and (r2,M*2,C'2.K'2)n8 using data from a set (Y',C',M\K')n2 where y>8 
and z<8 and where t>oth the "ny" and "nz" sets are the sets which produced a zero error closest to the set 

35 n = 8. 

The following illustrative example is provided to explain the optimization process based on a fixed 
number of value sets discussed above. 

Assume that the transform table 128 is to be constructed using the (YMCK)n sets of color values from 
Table II. Let us further assume that we are at a stage in the process where the sets of (YMCK)n shown in 

40 Table III have been sequentially introduced over line 154 to model 140, line 156 to adder 160 and line 155 
to LUT compiler 128. As a result, corresponding sets of (Y'M'C'K') of values and associated error En values 
have been obtained. The errors En at this time have not been tested for optimum error condition, for errors 
different than zero. The number of sets used for optimization will be limited to four (that is x = 4), to make 
the illustration simpler. The four sets of values and associated error values held In memory at the beginning 

45 of the optimization are shown In Table III. 



Table III 



50 



Set# 


(Y,M,C,K) 


(YMV1\C',K) 


E 


n 


76. 


76, 


76, 


89 


73, 


72. 


78, 


90 


0 


n + 1 


76. 


76. 


102, 


89 


73, 


73, 


100, 


89 


0.1 


n + 2 


76, 


76, 


128, 


89 


73. 


75, 


131, 


91 


0 


n + 3 


76, 


76, 


153, 


89 


73, 


71. 


152. 


93 


0 



55 

For input sets n, n + 2 and n + 3 the error is zero. Thus, we have perfect match and can accept the 
(Y.M.C.K) and corresponding (Y',M',C',K') values in the compilation of the transfer table 128. However, we 
need to check the results of set n + 1 to see if a better match may be available. 
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We proceed as we did earlier in the optimization process discussed for the (R.G.B) example. However, 
the process is done with two sets of input values to model 142, the sets that gave a zero match located on 
either side of the set that is optimized. In this example, the set that is optimized is set n + 1 and the two 
sets, that gave zero error on the two sides of set n + 1, are the n set and the n + 2 set. We apply the 

5 (Y.M.C.K)n+i (i.e., 76, 76, 102, 89) set of values again to model 140, but we apply to model 142 through the 
adder 160 first the (YM\/l',C',K')n (i.e., 73, 72. 78, 90) set of values which is the set corresponding to the "n" 
input, the set which gave a zero error in matching the Input values for "n". 

The (Lab)n+i and (L'a'b*)n outputs from the two models are compared in comparator 166 and new 
correction factors dY'n. dM'n, dC'n. are generated which result in a new set. Y'1n+i. M'1„+i, C'1n+i, K'1n+i of 

10 values when added to the (Y',M*,C',K*)n in adder 160. As before, a new set of Ua'b* values is produced and 
again compared in the comparator 166, and the process repeated, including the K incrementing step, until a 
minimum error value is found again. Assume these new values to be Y'1n+i=73. M'1n+i=74, C'1n+i =100 
and K'1 n+t —88. Let us further assume that this set of values results in a perfect match and a zero error. We 
temporarily store these results. 

15 Next we repeat the above process using this time the set of (Y*M'C'K*)n+2 values as an input to the 
adder 160, while the (Y,M,C,K)n+i set continues to be applied to model 140. This time we find that another 
set of (Y!2.M'2.C'2.K*2.)„+i values, i.e., r2n+i=71. M'2n+i=76. C'2„+i = 100 and K'2n+i=91, where the 
error is again zero. These values are also stored. These result are shown In Table IV: 

20 TABLE IV 



Set# 


(Y.M.CK) 




(YMVI 








E 


n 


76, 


76, 


76, 


89 




73. 


72, 


78, 


90 


0 


n + 1 


76. 


76, 


102, 


89 


(a) 


73, 


73, 


100, 


89 


0.1 












(b) 


73, 


74. 


100, 


88 


0 












(c) 


71, 


76, 


100, 


91 


0 


n+2 


76, 


76, 


128, 


89 




73, 


75, 


131. 


91 


0 


n+3 


76. 


76. 


153, 


89 




73. 


71. 


152, 


93 


0 



30 

The question now becomes which set of values (Y\M',C',K')n+i should we keep in compiling our table 
128 for the input set (Y,M,C,K)n+i. Obviously set (a) is the worst match, and we can discard it. Of the other 
two, it Is preferred to keep the set that will provide the most accurate interpolation results, when 
interpolation involving those values must be used. One way to do this is to select the set that has values 
best distributed between the preceding "n" set and the following "n+2" set A way to make this selection is 
as follows. Let Table V represent the two possibilities of selection. 

TABLE V 



Set# 


Possibility #1 (Y\M\C',K) 


Possibility #2 (Y'.M'.C'.K') 


Input (Y.M.CK) 


n 


73, 


72, 


78, 


90 


73, 


72, 


78. 


90 


76, 


76. 


76, 


89 


n + 1(b) 


73, 


74. 


100, 


90 










76. 


76. 


102. 


89 


n + 1(c) 










71, 


76, 


100. 


88 


76, 


76. 


102, 


89 


n + 2 


73, 


75. 


131, 


91 


73, 


75. 


131. 


91 


76, 


76, 


128, 


89 



We next compare the point spread deviation between input (Y.M.G.K) sets and the corresponding 
(Y*,M;,C*,K) sets. For instance, comparing the inputs between sets n and n + 1, we have: Yn+i-Vn = 0, 

50 Mn+i=Mn = 0, Cn+i-Cn = 26, Kn+rKn = 0 and comparing the inputs between sets n + 1 and n + 2, we have: 
Yn+2-Yn+i=0, Mn+2TMn+i=0, Cn+2-Cn+i = 26. Kn+2-Kn+i=0. Next. WO do the Same using the values for 
Y',M'.C\K' from the first selection column (labelled "Possibility #1), and we get: Y'n+i-Y'n = 0. M'„+i-M'n=2. 
C'n+i-C*n = 22, K'n+i-K'n = -1. Thus. WO See that for a zero change in the Y input values, there Is a zero 
change in the Y' values. For zero change in the M input values, we have a change of +2 points in the M* 

55 value. For a 26 point change in the C value input, the C value changes by 22 points, or a short fall of -4. 
And for a zero input change in the K value we have a -1 point change in the K'. 

Similarly the differences observed between the n + 2 and n + 1 sets, yields a difference for Y = 0 points, 
M= +1. C= +5and K= +2. 
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We will refer to these differences in the values as the deviation from the ideal situation where the 
increments of change between input values are the same as the increments of change between correspond- 
ing values, that Is where (Cn-Cn-i) = (C'„-C'n-i). The average of the deviation magnitude AP„ for the 
deviation between sets n and n + 1 for Y\M\C\K' from Table V selection #1 column is ADi =7/4 or 1.75; for 
5 sets n + 1 and n + 2, again for selection #1 the average is AD2 =2. 

For values from column selection #2 in Table V, the same process yields AD3=2.75 for sets n and 
n + 1, and 2.25 for AD* for sets n + 1 and n + 2. Thus, after comparing the average deviations for the two 
selections, set n + 1(b) In selection #1 is seen to yield set values that have magnitudes which, when taken 
as adjacent pairs of values with the preceding or following sets of values define color spans that on the 
10 average deviate less from the spans defined by the corresponding input color pairs than the values 
obtained from the sets of values from selection #2. Therefore, the set of (YM\/l\C\K')n+i(b) values is used in 
the transform table compilation. If set (a) could not be discarded as obviously the worst match, it could be 
compared to sets (b) and (c) like set (b) was compared to set (c) above to determine which set is the best 
match. 

75 Figure 9 shows a situation where image color matching is desired in images produced with a color 
printer 214 which accepts RGB values (through a converter 235 that converts RGB inputs to YMCA using a 
mathematical relationship that the printer manufacturer has determined as an accurate color conversion) 
and a CRT monitor 212. For this arrangement, as shown in Figure 9, the output of work station 10 is an 
RGB signal over line 216. It is used to drive the CRT display 212 directly and the color printer 214 through 

20 a 3 to 4 color conversion unit 225 which converts RGB inputs to YMCK outputs; an adaptor 222 comprising 
a transform LUT 228 constructed in accordance with this invention is placed ahead of conversion unit 225. 

This LUT 228 is compiled using a process similar to the process illustrated in Figure 8. However, in 
generating the model 44 for the printer 214, the converter unit 225 is included with the color printer 214 as 
an integral part thereof. Thus, this model 44 which is schematically shown in generic format in Figure 1 1 

25 has an input 45' which leads to converter 47. The output of converter 47 over line 45 is directed to the input 
of a receiving and determining element 46 and therefrom in a similar arrangement to that shown in the 
generic model of Figure 4 produces an output 52 either directly through LUT 48, or by. Interpolation through 
interpolator 50. 

In the arrangement shown in Figure 10, the output of the work station 10 is a YMCK signal which drives 

30 directly color proofer 312 over lines 316 and 318. The YMCK output is also used to display an image on a 
CRT monitor 314. This monitor 314 is able to accept YMCK input values because of a 4 to 3 color 
conversion that occurs in a color converter 325 prior to being used by the monitor 314. To obtain more 
accurate image color matching, an adaptor 322 which includes a receiving and determining element 326 
and a transform LUT 328 for converting YMCK input values to Y'M'C'K' Is placed ahead of the converter 

35 325. As before an interpolator 330 is Included In the adaptor to derive values not present in the LUT 328. 
The Y'M'C'K' values are directed to a CRT display 314 over line 324 through existing converter 325. 

The transform LUT 328 is also generated using the process described earlier, wherein two models 44 
are created representing the two display devices 312 and 314. The CRT display device model 314. 
however, as in the previous example. Includes the 4 to 3 converter 325 and Is also similar to the generic 

40 model illustrated In Figure 11. The LUT 328 Is again compiled using the process Illustrated In Figure 8. The 
transform LUT 328 converts input YMCK values to output Y'M*C'K* which when sent to converter 325 
generate RGB values such that when displayed on the CRT display 314 will produce an Image having the 
same or substantially the same color appearance as the image produced by the color proofer 312 using the 
original YMCK values directly. 

45 In discussing this invention there has been constant reference to "interpolation". Interpolation leading to 

derived output sets of color image values based on input sets of image values is used in performing the 
process of this Invention. It is evident that if interpolation is not used, the LUTs employed both in the 
models and in the adaptors must supply values for all possible combinations of input color values. In a 
digital system operating in an 8bit domain, for a three variable system. I.e., RGB, the possible input value 

50 combinations are 256*256*256 or 16,777,216 combinations. To develop the necessary LUTs, each of these 
combinations would have to be printed as a patch, in accordance with this process, measured, compared, 
and optimized, clearly an impractical process because of the shear magnitude of numbers and calculations 
involved. Thus, limited numbers of sets are selected as illustrated in Tables I and II and interpolation is 
used for intermediate sets of input color values. 

65 The particular interpolation process used may vary. Depending on the accuracy one wishes in the 
system more or less calculation intensive interpolations may be used. In practicing this invention, It is 
preferred that a less calculation intensive Interpolation method is used In the model Interpolators, while a 
more sophisticated method is used in the adaptor interpolators. Examples of both types are given below. 

12 
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However, any of the interpolation methods described herein below, as well as other interpolation formulae, 
such as quadratic interpolation, can be used for either the adaptor interpolators, or the model interpolators, 
the final choice being based on desired accuracy of the matching process, and the availability of resources. 
This invention Is not restricted to the use of one interpolator for adaptors 22, 122, 222, 322 and a different 

5 one for the models 40, 42, 140 and 142 even though such arrangement is both used and preferred. 

The preferred interpolation performed in interpolators 30, 130, 230 and 330 shown in the adaptors 
shown in Figures 1. 6. 9 and 10, respectively, that is the interpolator 30, 130. 230 or 330 that accompanies 
a transfonm LUT 28, 128. 228 or 328, respectively. Is tetrahedral Interpolation performed through a 
programmed computer. A pair of appropriate programs for programming such a computer to perform this 

10 interpolation for three and four variables are included with this application as Appendices "A" and "B" 
which appear immediately before the Claims. 

Linear Interpolation, tri-linear or quadralinear depending on the particular color space In which the color 
values exist. Is the preferred choice of Interpolation in models 40, 42, 140 and 142 shown in Figures 3 and 
8. used in deriving the transform LUTs 28, 128, 228, 328, as it is simpler, and often adequate to produce 

75 good results. The algorithm for this type of interpolation as applied in this invention is next explained for the 
more complex case of a four to three transform, i.e., from (YMCK) Input values to (Lab) as would be used in 
models 140 or 142. The three to three (RGB color space) and four to four (YMCK color space) transform 
cases are similar. 

The numerical values for Y.M,C,K given in Table II represent data points along a system of axes Y, M, 
20 C. and K. In the case of three variables, i.e., R,G,B color values, this may be an orthogonal axis system that 
may be easily visualized. NATith four variables visualization Is very difficult, however, such actual visualization 
is not necessary to perform the interpolation. 

The data points representing predetermined color values of a variable along the axis for this variable are 
separated by a specified number of units representing a distance measured along this axis. Thus, the data 
25 points Y = 102 and Y=128 are separated by a distance of 26 units along the Y axis, while the data points 
M = 128 and M = 153 are separated by 25 units along the M axis. An input value may fall either on a data 
point, or between data points. When it falls on a data point the LUT provides a corresponding value. If the 
input value fall between data points, then an output Is calculated as follows. 

Let F, be the ratio of the distance of the input value from a preceding data point along an axis to the 
30 distance between that data point and the next consecutive data point along this axis. Thus, If the color value 
set has four components, (i = 0, 1 , 2, or 3, for the Y, M, C, K components) and Yn is a data point along the Y 
axis, Yn+i is the next data point along this axis, and Yi is the Y component value of a color Input set of 
values whose numerical value lies between Yn and Y„+i, then: F| = (Y1-Yn)/(Y„+i-Yn) for 1 = 0. 
We also define a quantity Qi = 1-Fi, and the coefficients tj (where j = 1 .... 16) as follows: 

35 

ti = Fo*F, *F2*F3 te = Qo*Qi ^FsTg ti i = Fo*Fi *Q2*Q3 

t2 = Qo*Fi *F2*F3 t7 = Qo*Fi *Q2T3 ti 2 = Fo"Ql *Q2*Q3 

t3 = Fo*Qi T2T3 ts = Qo"Fi *F2*Q3 ti 3 = Qo*Fi *Q2*Q3 
t4 = Fo*F, *Q2*F3 t3 = Fo*Qi ^Qz^Fa ti ♦ = Qo*Qi ^2*03 
40 t5 = Fo*Fi *F2*Q3 ti 0 = Fo*Qi *F2*Q3 ti 5 = Qb*Qi ^CfeTa 
ti6 =Qo*QrQ2*Q3 

Let L(Yn.Mn.cn,Kn). a(Yr>.Mn.cn.Kn). b(Yn.Mn.cn.Kn). roprosont the Lab values from the LUT for the model 140 for 
instance, for an Input Y„.M„,Cn.Kn, L(Yn+i.Mn.cn.Kn). a(Yn+i .Mn.cn.Kn). b(Yn+ 1 .Mn.cn.Kn) tho Lab values for an Input 
45 Yn+1. Mn, Cn, and Kn, L(Yn.Mn+ 1 .cn.Kn), a(Yn.ivin+ i.cn.Kn). b(Yn.Mn+ 1 .cn.Kn). the values from the LUT for an input Yn, 
Mn+1. Cn and K^, and so on for all 16 combinations of input values (Y, M. C. K)n and (Y. M, C. K)n+i. The 
Lab values Li, ai, and bi for an input set of Y1,M1,C1,K1 color values falling between the n and n + 1 sets 
of values in the LUT. are given by the following equations: 

50 Li = ti * L(Yn+1.Mn+l,Cn+1.Kn+1) + t2 * L(Yn.Mn+1,Cn+1.Kn+1) 

+ t3 * L(Yn+1.Mn,Cn+1.Kn-i-1) + t4 * L(Yn+1,Mn+1,Cn,Kn+1) 

+ ts " L(Yn+1,Mn+1,Cn+l.Kn) + t6 * L(Yn.Mn,Cn+ 1.Kn+ 1) 

+ t7 * L(Yn,Mn+l.Cn.Kn+T) + t8 * L(Yn.Mn+ 1.Cn+ l.Kn) 

+ to * L(Yn+l,Mn,Cn,Kn+1) + ti 0 * L(Yn+ l.Mn,Cn+ l,Kn) 

55 + til * L(Yn+l.Mn+1,Cn,Kn) + ti 2 * L(Yn+ l.Mn,Cn.Kn) 

+ t13 * L(Yn,Mn+1.Cn,Kn) + t1 4 * L(Yn,Mn.Cn+ 1 .Kn) 

+ t15 * L(Yn,Mn,Cn.Kn+1) + t16 * L(Yn,Mn.Cn,Kn)- 
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Similarly, 

ai = ti * a(Yn+l.Mn+1,Cn+l.Kn+1) + t2 * a(Yn.Mn+1.Cn+1.Kn+l) 

+ t3 * a(Yn+1.Mn.Cn+l.Kn+1) + t4 * a(Yn+ l.Mn+ 1.Cn,Kn+ 1) 

+ t5 * a(Yn+1,Mn+1,Cn+1.Kn) ^ t6 * a(Yn.Mn.Cr>+ 1,Kn+ 1) 

+ t7 * a(Yn.Mn+1.Cn,Kn+1) + t8 " a(Yn,Mn+ 1 .Cn+ 1 .Kn) 

+ t9 * a(Yn+1.Mn,Cn,Kn+1) + t1 0 * a(Yn+ 1,Mn.Cn+ 1,Kn) 

+ t11 * a(Yn+1.Mn+1.Cn,Kn) "** t1 2 * a(Yn+ I.Mn.Cn.Kn) 

+ t13 * a(Yn,Mn+1,Cn,Kn) + t1 4 * a(Yn,Mn.Cn+ 1 ,Kn) 

+ t15 * a(Yn.Mn.Cn.Kn+1) + t16 * a(Yn,Mn,Cn.Kn)- 

And. 

bl = ti * b(Yn+1.Mn+1,Cn+1.Kn+1) + t2 * b(Yn,Mn+ 1 .Cn+ 1 .Kn+ 1) 

+ t3 * b(Yn+1,Mn,Cn+1.Kn+1) + t4 * b(Yn+ 1.Mn+ 1.Cn.Kn+ 1) 

+ t5 * b(Yn+1,Mn+1.Cn+1,Kn) + t6 * b(Yn,Mn.Cn+ 1,Kn+ 1) 

+ t7 * b(Yn.Mn+1,Cn,Kn+1) + t8 * b(Yn,Mn+ 1.Cn+ l.Kn) 

+ t9 * b(Yn+1.Mn.Cn.Kn+1) + t10 * b(Yn+1,Mn.Cn+ 1 ,Kn) 

+ t11 * b(Yn+1.Mn+1,Cn.Kn) + ti 2 * b(Yn+ 1 .Mn.Cn.Kn) 

+ t13 * b(Yn,Mn+1.Cn.Kn) + 4 * b(Yn.Mn.Cn+ l.Kn) 

+ t15 * b(Yn.Mn,Cn.Kn+1) + t16 * b(Yn.Mn.Cn.Kn). 

Those skilled in the art, having the benefit of the teachings of the present invention as herein above set 
forth, can effect numerous modifications thereto. These modifications are to be construed as being 
encompassed within the scope of the present invention as set forth in the appended claims. 
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APPENDIX A 



*****4r**:fr ************************************ 

(c) Copyright 1992 
£. I. DuPont de Nemours and Conpany (Inc.) 
Wilmington, Delaware 19898. 

t.r3us_t . c 

A function for tetrahedral type highspeed interpolation 
for performing 3 to 3 or 3 to n color coordinate 
conversions for the case of uniform sampling along 
coordinate axes. The vertices of the ci^e are 
numbered sequentially in a natural binary order 
(x3 x2 xl) . 

C. B. Chittineni 

June 1, 1992 



tr3us_t (n,nsrxl,x2,x3,l,x, y, z, val) 
no of outputs (int short) 

total number of data sample points per output (int short) 
data point with in the cube along xl axis (int short) 
data point with in the cube along x2 axis (int short) 
data point with in the cube along x3 axis (int short) 
length of cube (int short) 

data vector at current node (least rapidly varying) 
(unsigned char) 

data vector at current node (second least rapidly 
varying) (unsigned char) 

data vector at current node (most rapidly varying) 
(unsigned char) 
val- pointer to interpolated values (n) (int short) 

Modification: (Chittineni: July 9^ 1992) Changed the vertex numbering 
as (n) -> (X y z) (0 - 1 1 1, 1-011, 2-101, 
3-00 1, 4-11 0, 5-01 0, 6-10 0, 
7-000) (Liang's data files format) • 

*********it*****ik,**1tii*iii[iiii**it**1t*1t1r*1r *********** / 

tr3us_t(n, ns, xl, x2, x3, 1, x, y, z, new) 
short new [3], n, ns, xl, x2, x3, 1; 
unsigned char *x, *y, *z; 
{ 

int short wO, wl, w2, w3, a, b, c, d; 

int itm, j, jl, dt; 

int j2, j3, j4; 

unsigned char *x0, *yO, *zO; 

itm«0; 

if (x2 >- xl) { 

if (x3 >- x2) { 

itm-3; 
}else if (x3 >- xl) { 
itm-5; 
)else{ 

itm-4; 

} 

)else if (x2>-x3) { 



Module : 
Purpose : 



Author : 
Date Create: 



Usage : 

where n 

ns - 

xl - 

x2 - 

x3 - 

1 - 

X - 

y - 
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itm»0; 
}else if (x3 >- xl) { 
±tin-2; 
}else{ 

itm-1; 

} 

switch (itm) { 
case 0: 

wO-(l-xl); wl-(xl~x2); w2-(x2-x3); w3-x3; 
new[0]-(wO*x(7]+wl*xI3]+w2*x[lJ+w3*x(0) ) /I; 
new[l]-(wO*y t7]+wl*y [3]+w2*y [l]+w3*y(0] ) /I; 
new 12)- (wO*z [7J+wl*2 [31+w2*2 (l]+w3*2 [0] ) /I; 
break; 

case 1: 

wO-(l-xl); wl-(xl"x3); w2-(x3-x2); w3-x2; 
new[0]-(wO*x[7]+wl*x[3J+w2*xI2]+w3*x[0])/l; 
newtl]-(w0*y[7]+wl*y(3)+w2*y(2]+w3*y[0]) /I; 
new[2]-(wO*z[7]+wl*z(3]+w2*z[2]+w3*2(0)) /I; 
break; 

case 2 : 

w0-(l-x3); wl-(x3-xl); w2-(xl-x2); w3-x2; 
new [ 0 ) - ( wO *x [ 7 ] +wl *x [ 6 ] +w2 *x [ 2 ] +w3 *x [ 0 ] ) / 1 ; 
new [1] - {wO*y [7] +wl*y [6] +w2*y [ 2 ] +w3*y [ 0 ] ) /I; 
newt21-(w0*zt7]+wl*z[6]+w2*zt2]+w3*z[0] )/l; 
break; 

case 3 t 

w0«(l-x3); wl-{x3-x2); w2-(x2-xl); w3=xl; 
new[0]-{w0*x[7)+wl*x[6]+w2*x[4]+w3*x(0) ) /I; 
new [1] - (wO*y 17) +wl*y [6] +w2*y [4] +w3*y [0] ) /I; 
new[2]-(w0*z [7)+wl*z [6]+w2*z [4]+w3*z [0] ) /I; 
break; 

case 4: 

w0»(l-x2); wl«(x2-xl); w2-(xl-x3); w3=x3; 
new[0]-(w0*x[7]+wl*x[5]+w2*x[l]+w3*x[0] ) /I; 
new [ 1] - (wO*y [7] -fwl*y [5) +w2*y [ 1 ) +w3*y [ 0] ) /I; 
newi'2]-{w0*z[7]+wl*z [5) +w2*z [ 1] +w3*z [ 0] ) /I; 

break; 

case 5: 

w0=(l-x2); wl-(x2-x3); w2«(x3-xl); w3-xl; 
new {0] - (wO*x [7) +wl*x [5) +w2*x [ 4 ] +w3*x [0] ) /I; 
new(l]-<w0*y[7]+wl*y [5)+w2*y [4]+w3*y [0] ) /I; 
new[2]-(w0*z[7)+wl*z[5)+w2*z[4)+w3*z(0)) /I; 
break; 

} 

return; 

} 
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APPENDIX B 



/*«*-*^*«**«*****iri«r* ************** ******* ************** 

(c) Copyright 1992 
£. I. DuPont de Hemours and Company (Inc.) 
Wilmington, Delaware 19898. 



Module: 
Purpose: 



Author : 



. tr4us_t.c 

A function for tetrahedral type highspeed interpolation 
for performing 4 to 4 or 4 to n color coordinate 
conversions for the case of uniform sampling along 
coordinate axes. The vertices of the hypercvtbe are 
numbered sequentially in a binary order (xl x2 x3 x4). 
It is to be noted that the ordering is in^ortant 
in arranging the saxRpled data and in passing the 
arguments to the function. 

C. B. Chittineni 



Date Create: June 1, 1992 

Usage: tr4us_t (n^ns^ xl,x2#x3f x4,l/ cblrmblr ybl, Icbl, val) 

where n - no of outputs (int short) 

ns - total number of data Scunple points per output (int) 
xl - data point with in the hypercube along xl axis (int short) 
x2 - data point with in the hypercube along x2 axis (int short) 
x3 - data point with in the hypercube along x3 axis (int short) 
x4 - data point with in the hypercube along x4 axis (int short) 
1 - length of hypercube (int short) 

cbl- data vector at current node (least rapidly varying) 
(unsigned char) 

mbl- data vector at current node (second least rapidly varying) 
(unsigned char) 

ybl- data vector at current node (third least rapidly varying) 

(unsigned char) 
kbl- data vector at current node (most rapidly varying) 

(unsigned char) 
val- pointer to interpolated values (n) (int short) 

Modification: (Chittineni : July 9, 1992) Changed the vertex niombering 
as (n) -> (cbl mbl ybl kbl) (0-1111, 1-0111, 
2-1011, 3-110 1, 4-1110, 5-0011, 
6-010 1, 7-0110, 8-100 1, 9-1010, 
10- 1 1 0 0,11 - 1 0 0 0,12 - 0 1 0 0,13 - 0 0 1 0, 
14- 0 0 0 1,15 - 0 0 0 0) (Liang's data format) . 

******* ************************************************* *************/ 

♦include <stdio.h> 
♦include <math.h> 

tr4us_t(n, ns, xl, x2, x3, x4, 1, cbl, mbl, ybl, kbl, new) 
int short n, xl, x2, x3, x4, 1; 
int short *new; 
int ns; 

unsigned char *cbl, *xribl, *ybl, ♦kbl; 

{ 



int short wO, wl, w2, w3, w4, a, b, c, d, e; 
int itm, j, jl, j2/ 



/* Locate the tetrahedran */ 
if (xl >- x2) { 
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if {x2 >- x3) ( 

i£(x3 >- x4) { 

itm-l; 
}else If (xl >- x4) { 
5 if (x2 >- x4) { 

itin-2; 
)else{ 

±tin-3; 

) 

}else ( 

10 itm-l 9; 

) 

)else if (x3 >- x4) { 

if (xl>-x4) { 

if (x2>-x4) { 

if (xl>«x3) { 

75 itin-5; 

)else{ 

itin-13; 

) 

}else if (xl >- x3) { 
itin-6; 
}else{ 

20 itm-l 4; 

} 

}else{ 

itm-15; 

) 

}else if (xl >- x4) { 
25 itro-4; 

)else if (xl >- x3) { 

itm-20; 
}else { 

itm-21; 

} 

}else if (x2 >- x3) { 

if (x3 >- x4) { 

if (xl >- x4) { 

if (xl >- x3) { 

itm-7/ 
}else { 

itm-l 1; 

35 ) 

}else{ 

itm-12; 

) 

}else if (xl >- x4) ( 
itm-8; 

40 )else if (x2 >- x4) { 

if(xl >« x3){ 

itm-9; 
}else{ 

itm-10; 

} 

45 }else if(xl >« x3) { 

itm'»23; 
}else{ 

itm=24; 

}else if (x3 >- x4) { ^ 
if (xl >- x4) { 

50 itm-l 7; 

)else if (x2 >- x4) { 
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ltm-18; 
}else{ 

itin-16; 

} 

>else{ 

itin-22; 

} 

fprintf (stderr, "%d \n", itm);*/ 
switch (itm) { 
case 1: 

wO-(l-xl); wl-{xl-x2); w2-(x2-x3); w3-(x3-x4); w4-x4; 

new 1 0 ] - (wO*cbl [15] +wl*cbl 111) +w2*cbl [ 10 ] +w3*cbl [ 4 ] +w4*cbl 1 0 ] ) /I 

new [ 1 ) - (wO*mbl [ 15 ) +wl*xnbl { 11 ) 4-w2*inbl [ 10 ) +w3*inbl [ 4 ) +w4 *inbl 1 0 ] ) /I 

new [2 ] - (wO*ybl [ 15) -i-wl*ybl [ 11 ) +w2*ybl ( 10 ) +w3*ybl [ 4 ) +w4*ybl [ 0 ) ) /I 

new 1 3 ] - ( wO *)cbl [15] +wl*kbl ( 1 1 ) +w2 *kbl [ 10 ) +w3*)cbl 1 4 ) +w4 *)cbl [ 0 ] ) / 1 

brea]c; 

case 2: 

wO-(l-xl); wl-(xl-x2); w2-(x2-x4); w3-(x4-x3); w4-x3; 
new[0]-(wO*cbl[15)+wl*cbl{ll)+w2*cbX(10)+w3*cbl[3]+w4*cbl[0] ) /I. 
- new[l]-(wO*mbl[15)+wl*xhblIll)+w2*ntolI10)+w3*itiblI3)+w4*inbl[0] ) /I. 
new[2]-(w0*ybl[15)+wl*ybl[ll)+w2*ybl[10]+w3*ybl[3]+w4*ybl[0] ) /I. 
new[3]-(w0*]cbl[15]+wl*)cbl[ll)+w2*kbl[10)+w3*kbl[3]+w4*3cbl[0) ) /I. 
brealc; 

case 3: 

wO-(l-xl); wl-(xl-x4); w2-(x4-x2); w3"(x2-x3); w4-x3; 

new [ 0 ] - (wO*cbl ( 15 ) +wl*ebl [ 11 ] +w2*cbl [ 8 ] +w3*cbl [ 3 ) +w4*cbl [ 0 ) ) / 1 ; 

new [ 1 ] - (wO*mbl ( 15) +wl*mbl [ 11 ) +w2*inbl ( 8 ) +w3*iribl [ 3 ) +w4 *itibl C 0 ] ) /I ; 

new[2]-(w0*ybl(15)+wl*ybl[ll)+w2*ybl[8)+w3*yblC3]+w4*ybl[0] ) /!;. 

ne w [ 3 ] - ( wO*)cbl [15] +wl*kbl [11) +w2 *)cbl [ 8 ) +w3 *)cbl [ 3 ] +w4 *3cbl ( 0 ) ) / 1 ; 

brea)c; 

case 4: 

wO-(l-xl); wl-(xl-x4); w2-(x4-x3); w3-(x3-x2); w4-x2; 
new(0]-(w0*cbl[15)+wl*cbl[ll)+w2*cbl[8]+w3*cbl[2)+w4*cbl[0] ) /I; 
ne w [ 1 ] - < wO *inbl [ 1 5 ] +wl *nibl [ 1 1 ) +w2 *inbl [ 8 ] +w3 *rnbl [ 2 ) +w4 *inbl [ 0 ] ) / 1 ; 
new 1 2 ] - < wO *ybl [15] -i-wl *ybl [11] +w2 *ybl [ 8 ] +w3 *ybl [ 2 ) +w4 *ybl [ 0 ] ) / 1 ; 
new ( 3^ - (wO*]cbl [ 15 ) +wl*Jcbl [ 11 ) +w2*)cbl [ 8 ) +w3*)cbl [2 ) +w4 *kbl 10}) /I; 
breaks- 
case 5: 

wO-(l-xl); wl-(xl-x3); w2"-{x3-x2); w3-(x2-x4); w4-x4; 

new [ 0 ] - ( wO*cbl [ 15 ) +wl*cbl [11] +w2*cbl ( 9 ] +w3 *cbl [ 4 ] +w4 *cbl [ 0 ] ) /I ; 

new [ 1 ] - (wO*rnbl [ 15] +wl*robl [ 11 ) +w2*inbl [ 9] +w3*inbl [ 4 ] +w4 *inbl ( 0 ] ) /I; 

new ( 2 ] - ( wO*ybl [ 15 ] +wl*ybl [ 11 ] +w2*ybl ( 9 ) +w3*ybl [ 4 ] +w4 *ybl [ 0 ) ) / 1 ; 

ne w [ 3 ] - ( wO*kbl (15) +wl*kbl [ 1 1 ) +w2 *kbl [ 9 ) +w3*]cbl [ 4 ] +w4 *kbl [ 0 ) ) / 1 ; 

break; 

case 6 : 

wO=(l-xl); wl-(xl-x3); w2-(x3-x4); w3-(x4-x2); w4-x2; 
new[0]-(w0*cbl[15)+wl*cbl[ll]+w2*cbl[9]+w3*cbl [2] +w4*cbl [0] ) /I; 
new[l]-(w0*inbl(15]+wl*inbl[ll]+w2*inbl(9]+w3*mbl [2] +w4*inbl [0] ) /I; 
new ( 2 ] - (wO*ybl [ 15] +wl*ybl [ 11 ] +w2*ybl [ 9] +w3*ybl (2 ] +w4 *ybl [ 0 ] ) /I; 
new ( 3 ) = ( w0*kbl [15) +wl *kbl [11] +w2 *kbl [ 9 ] +w3 *kbl [ 2 ] +w4 *kbl [ 0 ] ) / 1 ; 
break; 

case 7: 

w0-<l-x2); wl-(x2-xl); w2-(xl-x3); w3-(x3-x4); w4-x4; 
newI0]-(w0*cbl(15]+wl*cbl[12J+w2*cbl (10]+w3*cbl[4] +w4*cbl(0] ) /I; 
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new (1 ] - (wO*inbl [ 15] +wl*inbl 112] +w2*inbl 1 10 ] +w3*mbl 1 4 ) +w4 *inbl ( 0 ] ) /I ; 
new [2 ] - (wO*ybl C 15] +wl*ybl [ 12 ] +w2*ybl 1 10 ] +w3*ybl ( 4 ] +w4 *ybl ( 0 ) ) /I ; 
new [3 ] - (wO*)cbl ( 15] +wl*3cbl 1 12 ] +w2 *)cbl 1 10 ) +w3*kbl [ 4 ] -f w4*kbl [ 0 ] ) /I ; 
breaks- 
case 8: 

w0-(l-x2); wl-(x2-xl); w2->(xl-x4); w3-(x4-x3); w4-x3; 

new [0] - (wO*cbl [15] +wl*cbl [12)+w2*cbl I10]+w3*cbl 13] +w4*cbl [0] ) /I; 

new [ 1] - (wO*inbl 1 15] +wl*inbl 112] +w2*xribl ( 10 J+w3*mbl [3] +w4*inbl [ 

new [ 2 ] - (wO*ybl (15] +wl*ybl 1 12 ] +w2*ybl [ 10] +w3*ybl ( 3 ] +w4 *ybl ( 0 ] ) /I; 

new ( 3 ] - ( wO *3cbl [ 1 5 ] +wl *]cbl ( 1 2 J +w2 *Jcbl ( 10 ] +w3*kbl [ 3 J +w4 *)cbl ( 0 ] > /I ; 

break; 

case 9 : 

w0-(l-x2); wl-<x2-x4); w2-{x4-xl); w3-(xl-x3); w4-x3; 
new(0]-<w0*cbl[15]+wl*cbl(12)+w2*cbl[6]+w3*cbll3]+w4*cbl{0] ) /I; 
ne w [ 1 ] - { wO *inbl [15] +wl *inbl [ 12 ) +w2 *inbl ( 6 ] ■i-w3 *inbl ( 3 ) +w4 *inbl ( 0 ] ) / 1 ; 
ne w [ 2 ] - < wO *ybl [15] +wl *ybl [ 12 ] +w2 *ybl [ 6 ] +w3 *ybl ( 3 ] +w4 *ybl [ 0 ] ) / 1 ; 
new [ 3 ] - (wO*kbl [ 15] +wl*kbl [ 12 ] +w2 *kbl ( 6 ] +w3*kbl [ 3 ] +w4 *kbl [ 0 ]) /In- 
break; 

case 10 : 

w0-(l-x2); wl'-(x2-x4); w2-(x4-x3); w3«(x3-xl); w4-xl; 

new [03- ( w0*cbl [ 15 J +wl*cbl ( 12 ] +w2*cbl [ 6 ] +w3 *cbl [ 1 ] +w4 *cbl [ 0 ] ) /I; 

new [ 1 ] - ( wO *inbl [15] +wl *inbl [12] +w2 *iiibl [ 6 ] +w3 *mbl [ 1 ] +w4 *mbl [ 0 ) ) / 1 ; 

new[2]-(w0*ybl[15]+wl*ybl[12]+w2*ybl(6]+w3*ybl(l]+w4*ybl[0])/l; 

new [ 3 ] - ( wO *kbl [ 15 ] 4.W1 *kbl [ 12 ] +w2 *kbl ( 6 ] +w3 *kbl [ 1 ] +w4 *kbl [ 0 ] ) /I; 

break; 

case 11 : 

w0-(l-x2); wl-(x2«x3); w2-{x3-xl); w3-(xl-x4); w4-x4; 
new[0)-(w0*cbl[15]+wl*cbl[12]+w2*cbl[7]+w3*cbl [4]+w4*cbl [0] ) /I; 
new [ 1 ] - { wO *inbl [15] +wl *nibl ( 12 ] +w2 *iribl [ 7 ] +w3 *iribl [ 4 ] +w4 *inbl ( 0 ] ) / 1 ; 
new { 2 ] « ( wO*ybl [15] +wl*ybl [ 12 ] +w2 *ybl [ 7 ] +w3 *ybl [ 4 J +w4 *ybl [ 0 ] ) / 1 ; 
new ( 3 ] - ( wO *kbl [ 15 ] +wl*kbl [ 12 ] +w2 *kbl [ 7 ] +w3 *kbl ( 4 ] +w4 *kbl [ 0 ] ) /I; 
break; 

case 12 : 

w0-(l-x2); wl-(x2-x3); w2-(x3-x4); w3-(x4-xl); w4-xl; 
new[0]-(w0*cbl[15]+wl*cbl[12]+w2*cbl (7]+w3*cbl [1]-Kw4*cbl [0] ) /I; 
new [ 1 ] - {wO*inbl [15] +wl*mbl [ 12 ] +w2 *nibl [ 7 J +w3*nibl [ 1 ] +w4 *robl [ 0 ] ) /I; 
new ( 24 = <wO*ybl [15] +wl*ybl [ 12 ] +w2 *ybl ( 7 ] +w3*ybl ( 1 ] +w4 *ybl [ 0 ] ) /I; 
n ew [ 3 ] « { wO *kbl [ 1 5 ] +wl *kbl [ 1 2 ] +w2 *kbl ( 7 ] +w3 * kbl [ 1 ] +w4 * kbl [ 0 ] ) / 1 ; 
break; 

case 13: 

w0-(l-x3); wl-(x3-xl); w2-(xl-x2); w3-(x2-x4); w4-x4; 
new[0]-{w0*cbl[15]+wl*cbl[13]+w2*cbl[9]+w3*cbl[4]4-w4*cbl[0] )/l; 
new[l]*(w0*na5l[15]+wl*na3l[13]+w2*inbl[9]+w3*mbl[4]+w4*inbl[0] ) /I; 
new[2]-(w0*ybl[15]+wl*ybl[13]+w2*ybl[9]+w3*ybl[4]+w4*ybl[0] ) /I; 
new[3]-{w0*kbl[153+wl*kbl[13]+w2*kbl[9]+w3*kbl[4]+w4*kbl[0] ) /I; 
break; 

case 14 : 

w0-(l-x3); wl-(x3-xl); w2-(xl-x4); w3-(x4-x2); w4-x2; 
new[0]-(w0*cbl[15]+wl*cbl(131+w2*cbl[9]+w3*cbl[2]+w4*cbl [0] )/l; 
new[l]-<wO*^a^l[15]+wl*^a)l[13]+w2*^0bl[9] +w3*nibl I2)+w4*mbl [0] ) /I; 
new [ 2 ] - ( w0*ybl [ 15 ] +wl*ybl [ 13 ] +w2 *ybl ( 9 ] +w3 *ybl [ 2 ) +w4 *ybl [ 0 ] ) / 1 ; 

new[3]-(wO*kbl[15J+wl*kbl(13]+w2*kbl[93+w3*kbl[2)+w4*kbl[0] ) /I; 
break; 

case 15: 

w0-(l-x3); wl-<x3-x4); w2-<x4-xl); w3«(xl-x2); w4-x2; 
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new 1 0 ] - ( wO*cbl ( 15 ) +wl*cbl [ 13 ] +w2 *cbl ( 5 ] +w3 *cbl ( 2 ] +w4 *cbl ( 0 ] ) /I; 
new [ 1 ] - ( wO *inbl [15] +wl*rnbl [13] +w2 *inbl [ 5 ] +w3 *inbl ( 2 ) +w4 *mbl [ 0 J ) / 1 ; 
new[2]-(w0*ybl[15]+wl*yblI13]+w2*ybl[5]+w3*ybl(2]+w4*ybl(03)/l; 
newt3]-(w0*kbltl5]+wl*)cbl(13]+w2*kblI5]+w3*kbl[2]+w4*)cbl(03)/l; 

breaks- 
case 16: 

w0-{l-x3); wl-(x3-x4); w2-(x4-x2); w3-(x2-xl); w4-xl; 

new [ 0 ] - ( wO*cbl [ 15 ] +wl*cbl ( 13 ] +w2*cbl ( 5 1 +w3 *cbl [ 1 ] +w4 *cbl [ 0 ] ) /I; 

new [ 1 ] - (wO*inbl 1 15] +wl*robl [ 13] +w2*inbl [5] +w3*robl [1] +w4 *mbl 1 0] ) /I; 

new ( 2 ] - { wO*ybl [15] +wl*ybl 1 13 ] +w2 *ybl [ 5 ] +w3 ♦ybl { 1 ] +w4 *ybl 1 0 ] ) /I ; 

ne w [ 3 ] - ( wO*kbl 115] +wl*kbl 1 13 ] +w2 *kbl [ 5 1 +w3 *kbl ( 1 ] +w4 *kbl 1 0 1 ) /I ; 

break; 

case 17: . „ 

w0-(l-x3); wl-(x3-x2); w2-(x2-xl); w3-(xl-'x4); w4-x4; 
new [ 0 ] - ( wO*cbl { 15 ] +wl*cbl ( 13] +w2*cbl [ 7 ) +w3*cbl [4 ] +w4 *cbl 10]) /I; 
new 1 1 ] - (wO*iribl 1 15 ] -f-wl^inbl [ 13 ] +w2*inbl ( 7 ] +w3*inbl [4 ] +w4 *inbl 1 0 ] ) /I; 
new [ 2 ] - ( wO*ybl [15] +wl*ybl ( 13 ] +w2 *ybl [ 7 ] +w3*ybl ( 4 ] +w4 *ybl 1 0 ] ) /I; 
new [ 3 ] - ( wO *)cbl [ 15 ] +wl*kbl 1 13 ] +w2 *kbl ( 7 ] +w3 *kbl [ 4 ] +w4 *Wbl [ 0 ] ) /I ; 
break; 

case 18 : 

w0-(l-x3); wl-(x3-x2); w2-{x2-x4); w3-<x4-xl); w4-xl; 
new[01-(w0*cbl[15]+wl*cbl[13]+w2*cbl[7]+w3*cbl(l]+w4*cbl[0])/l; 
new [ 1 ] - ( w0*inbl 1 15 J +wl*inbl ( 13 ] +w2 *nibl [ 7 ) +w3 *robl [ 1 ] +w4 *nObl 1 0) ) / 1 ; 
new [ 2 ] - ( wO*ybl [ 15 J +wl*ybl 1 13] +w2*ybl [ 7 ] +w3*ybl [ 1 ] +w4 *ybl [ 0 ] ) / 1; 
newI3]-(w0*kbl[15)+wl*kblll3]+w2*kbl[7]-»-w3*kbl[l]+w4*kbl(0])/l; 

break; 
case 19 • 

w0-(l-x4); wl-(x4-xl); w2-(xl-x2); w3-(x2-x3); w4-x3; 
newI0]-(w0*cbltl5]+wl*cblI14]+w2*cblt8]+w3*cbl[3]+w4*cbl[0])/l; 
new ( 1 ] - ( wO *na>l [ 1 5 ] +wl*xnbl [ 14 ] +w2 *inbl C 8 ] +w3 *ihbl ( 3 ] +w4 *inbl 1 0 ] ) / 1; 
new[2]-(w0*ybl[15]+wl*ybl{14]+w2*yblt8]+w3*ybl[3]+w4*ybl[0]) /I; 
new [ 3 ] - (w0*kbl [ 15] -fwl*kbl [ 14 ] +w2*kbl [ 8 ] +w3*kbl f 3 ] +w4*kbl ( 0 ] ) /I; 
break; 

case 20: 

wO-(l-x4); wl-(x4-xl); w2-(xl-x3); w3-(x3-'x2); w4-x2; 

new [ 0 ] - ( wO *cbl [ 15 ] +wl*cbl [ 14 ] +w2 *cbl ( 8 ] +w3*cbl ( 2 ] +w4 *cbl ( 0 ] ) /I 

new [ 14^- ( wO *mbl (15) +wl *iTibl [14] +w2 *inbl [ 8 ] +w3 *inbl [ 2 ] + w4 *mbl [ 0 ] ) / 1 

new[2]-(w0*ybl[15]+wl*ybl[14]+w2*ybl[8]+w3*ybl(2]+w4*ybl[0] ) /I 

new ( 3 3 - ( wO *kbl ( 15 ] +wl*kbl [14] +w2 *kbl 1 8 ] +w3 *kbl 1 2 ) +w4 *kbl [ 0 ] ) / 1 

break; 

case 21 : 

w0-a-x4); wl-(x4-x3); w2-(x3-xl); w3"(xl-x2); w4-x2; 
new[0]-(w0*cbl[15]+wl*cbl[14]+w2*cbl(5]+w3*cbl(2]+w4*cbl[0] ) /I; 
new[l]-<w0*xnbl[15)+wl*nvbl[14]+w2*xnblI5]+w3*mbl[2)+w4*mbl[0] ) /I; 
new[2]-(w0*ybl[15]+wl*ybl[14]+w2*ybl[5]+w3*ybl[2]+w4*ybl[0] )/l; 
new [ 3] - <wO*kbl ( 15] +wl*kbl [ 14 ] +w2*kbl [5] +w3*kbl [2 ] +w4 *kbl ( 0 ) ) /I; 
break; 

case 22: 

w0-(l-x4); wl-(x4-x3); w2-(x3-x2); w3«(x2-xl); w4-xl; 
new(0]-(w0*cbl[15]+wl*cbl(14]+w2*cbll5]+w3*cbl[l]+w4*cbl[0] ) /I 
new [ 1 ] - ( w0*inbl ( 15] +wl*mbl [ 14 ] +w2*mbl [ 5 ] +w3*rnbl [ 1 ] +w4*mbl [ 0 ] ) /I 
new[2) = (w0*ybl (15]+wl*ybl [14]+w2*ybl [5] +w3*ybl [1] +w4*ybl [0] ) /I 
new{3]-<w0*kbl[15]+wl*kbl(14]+w2*kbl[5]+w3*kbl(l]+w4*kbll0] ) /I 
break; 

case 23: 
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w0-(l-x4); wl-(x4-x2); w2-(x2-xl); w3-(xl-x3)/ w4-x3; 

new [ 0 1 - ( wO * cbl [15] -f wl * cbl [14] +w2 * cbl [ 6 ] +w3 * cbl [ 3 ] +w4 * cbl [ 0 ) ) / 1 

new [ 1 ] - ( wO *mbl [15] +wl *inbl [14] +w2 *inbl [ 6 ] +w3 *inbl ( 3 ] +w4 *inbl ( 0 ] ) / 1 

new(2]-(w0*ybl(15]+wl*ybl(14]+w2*ybl[6]+w3*ybl[31+w4*ybl[0]) /I 

new[3]-(wO*)cbl(151+wl*3cbl[14]+w2*3cbl[6]+w3*)cbl(3]+w4*kbl[0])/l. 

break; 

case 24: 

wO-(l-x4); wl«(x4-x2); w2-(x2-x3); w3-(x3-xl); w4-xl; 
new[0]-(w0*cbl[15]+wl*cbl(14]+w2*cbl[6]+w3*cbl[l]+w4*cbl [0] ) /I; 
new [ 1 ] - ( wO*iia3l [15] +wl*xnbl [ 14 ] +w2*inbl [ 6 J +w3*mbl ( 1 ] +w4 *nibl [ 0 ] ) /I ; 
n e w [ 2 ] - ( wO *ybl [15] +wl *ybl [14] +w2 *ybl ( 6 ] +w3 *ybl [ 1 ] +w4 *ybl [ 0 ] ) / 1 ; 
new[3]-(w0*)cbltl5]+wl*)cbl[14]+w2*Jcbl(6]+w3*kbl[l]+w4*Jcbl[0]) /I; 
break; 

) 

return; 

1 



Claims 

1. A method for matching the color display of at least a first device and at least a second display device, 
the first and second display devices having an input, comprising: 

(I) creating a transform LUT for converting input color values to output color values by: 

(1) producing a first preselected plurality of input color values; 

(2) using said plurality of input values to display a plurality of color patches in said first and said 
second display devices; 

(3) obtaining a colorimetric value of each of the displayed color patches in each of the displayed 
devices and using said colorimetric values to create a first and a second model for the first and 
second devices respectively correlating preselected input color values to displayed colorimetric 
values for each of the two devices; 

(4) inputting to said first and second models a second preselected plurality of color values; 

(5) comparing the output of the first model to the output of the second model to obtain an error 
signal indicative of the difference between the two output signals; 

(6) using the error signal to modify the input color values to the second model and again 
comparing the output of the first model to the output of the second model to obtain a new error 
signal; 

(7) repeating the process of steps (5) and (6) above until the error signal is a minimum; and 

(8) using the modified color values to create a transform LUT correlating input values to modified 
values; 

(II) using the transform LUT to transform any plurality of source color image values to corresponding 
new image values before inputting said new image values to the input of the second display device; 
and 

(III) displaying said new color image values on said second display device. 

2. The method according to Claim 1, wherein in performing the step II, first testing every one of the 
plurality of source color image values; and 

(a) if the source color image values are included in the LUT obtaining from the LUT new image 
values corresponding to the source color image values; and 

(b) If the source color image values are not included in the LUT obtaining by interpolation from 
adjacent source image values in the LUT new image values corresponding to the source image 

values. 

3. The method according to Claim 1 or 2. wherein between steps (l)(7) and (l)(8) there is the additional 
step of optimizing the minimum error. 



22 



nRn47fi.sAi I > 



EP 0 604 755 A1 



4. The method according to Claim 1 or 2, wherein prior to step (l)(5). the gamut of the output color values 
of the first model Is adjusted to be commensurate with the color values gamut of the second model. 

6. The method according to Claim 2, wherein the source color Image values are RGB values and the 
display devices are CRT displays. 

6. The method according to Claim 2, wherein the source color image values are YMCK values and the 
display devices are 4-color printing devices. 

7. The method according to Claim 2, wherein one of the display devices is a CRT display and one a 4- 
color proof er and wherein there is Included In the color proofer a color converter for converting RGB 
color values to YMCK color values. 

8- The method according to Claim 2, wherein one of the display devices is a CRT display and one a 4- 
color printer and wherein there is included in the CRT a color converter for converting YMCK color 
values to RGB color values. 

9. Apparatus for generating a transform LUT for converting Input color values to output color values. 

comprising: 

(1) digital color image values input means; 

(2) first and second display model means for producing each an output colorimetric tristimulus 
Image values for input color values, each of said model means having an input and an output, the 
second model Input connected to the input means; 

(3) an adder device having a first signal input connected to the input means, a second, correction 
error signal Input, a control signal Input, and an output connected to the first model input, for 
outputting modified color values; 

(4) means, connected to the output of said first and second model means, for comparing the 
tristimulus values output of the first and second display models and for producing an error signal; 

(5) means for testing the error signal to determine if said signal Is a minimum error signal, and for 
outputting a correction error signal and a control signal to the adder, said testing means connected 
between said means to compare and said adder device; and 

(6) means also connected to the input means, the adder output, and the means for testing, for 
receiving the adder output and for generating a transform LUT correlating digital input color values to 
the modified color values appearing at the output of the adder when the means for testing 
determines that the error signal is a minimum. 

10. The apparatus according to Claim 9. further comprising means for storing the generated LUT. 

11. The apparatus according to Claim 9, further comprising a gamut mapping means for mapping output 
values from the second model into a color gamut commensurate to that of the first model, the gamut 
mapping means positioned between the output of the second model and the means to compare. 

12. The apparatus according to Claim 9. further comprising an optimization means in the error testing 
means for determining if the minimum error is an optimum minimum error. 
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